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Abstract 



The effects of spectator light quark in decays of heavy mesons are considered, 
^ I which vanish in the hmit of factorization of matrix elements of four-quark operators 

over the mesons. These effects include the difference of the total widths as well as of 
the semileptonic decay rates between the and Ds mesons and also a contribution 
to inclusive semileptonic decay rates of and into the channel i v related 
to determination of the weak mixing parameter Vuh- If the observed difference of the 
lifetimes between Dg and mesons is attributed to non-factorizable terms, such terms 
can naturally give rise to a significant difference in inclusive semileptonic decay rates 
of these mesons, and to a light-flavor dependent contribution to decays B Xu^f- 
The latter contribution affects mostly the upper part of the inclusive spectrum of the 
invariant mass of the lepton pair, and may significantly exceed the previously claimed 
in the literature theoretical uncertainty in determination of from that part of the 

spectrum. 



1 Introduction 



The well known difficulty of determination of the weak mixing parameter \Vub\ from the 
inclusive decay B — > X„ £ z/ is in discriminating this process from the much more abundant 
decays B Xciiy (see e.g. in the review in Ref. [Q and references therein). The experi- 
mental method used so far is to impose a lower cutoff on the charged lepton energy at or 
above the kinematical boundary for the processes with charmed hadrons in the final state. 
Theoretically however a calculation of the lepton energy spectrum in the endpoint region is 
prone to a substantial uncertainty which is unlikely to be resolved without knowledge of 
B meson structure functions 0]. 

A somewhat more theoretically tractable scheme arises if a cutoff is applied to the 
invariant mass squared of the lepton pair, > q^, instead of the single- lepton energy. 
Clearly, the condition Qq > {Mb — M^)^ makes a sample of semileptonic decay events free 
from the background from B Xcd-v. The total decay rate into the part of the phase space 
limited by such constraint can be analysed, including nonperturbative contributions as well 
as the perturbative one, by applying the operator product expansion (OPE) for inclusive 
decay rates @, |^. Most recently an analysis along these lines was performed in Ref. 

The difficulty arising in this approach is that the "short distance" parameter in the OPE 
for the constrained inclusive decay rate is determined by the typical momentum /i^ of the 
u quark in the constrained kinematical region, where its maximal value is (M^ — go)/2-^s, 
which becomes slightly smaller than the mass of the charmed quark mc, even if the cutoff pa- 
rameter is chosen at exactly the kinematical limit for charmed final states, > {Mb — MoY- 
This puts the theoretical status of the OPE for the constrained inclusive rate on par with or 
even somewhat worse than that of an OPE based analysis of total inclusive weak decay rates 
of charmed mesons and baryons. It is well known that in the latter analysis the nonperturba- 
tive terms in OPE, whose relative contribution is suppressed as and m~^, are essentially 
as significant as the "leading" perturbative (parton) term, which behavior prominently man- 
ifests itself in the large differences of the lifetimes both between the charmed mesons and 
baryons and within the meson and baryon light-flavor multiplets. This simple observation 
illustrates the acute necessity of considering the higher terms, formally suppressed by ju~^ 
and in the OPE based analysis of the constrained decay rate oi B ^ Xud-V- The 
leading perturbative term as well as the nonperturbative contribution proportional to yU~^ 
were considered in Ref. , while the contribution relatively suppressed by was relegated 
to a theoretical uncertainty. The magnitude of the uncertainty, arising from this source. 
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was estimated (essentially by guessing) to correspond to about 5% theoretical error in \ Vub\- 
The purpose of this paper is to discuss the contribution of the third term in OPE which 
arises from non-factorizable terms in matrix elements of four-quark operators over heavy 
mesons, and to analyse, to an extent, the numerical significance of this contribution. The 
impact of non-factorizable terms on semileptonic decays of D and B mesons as well as on 
the difference of lifetimes of Dg and has been discussed previously 0, ||, § • It is however 
worthwhile, in view of the recent discussion P] of theoretical uncertainty in determination 
of \Vub\, to point out the effects of the non-factorizable terms in this particular issue, and 
to present specific estimates relevant to probing these terms from data on charmed mesons. 
It is emphasized here that these terms, generally different for and B"^ mesons (as first 
pointed out in Ref . ) , can be quite essential in the constrained inclusive semileptonic decay 
rates due to the h u transition, if the magnitude of these effects is estimated, following 
the arguments of Ref. [Q, from the observed difference of lifetimes of the Dg and mesons. 
A more direct evaluation of the light-flavor dependent part of the relevant non-factorizable 
terms would be possible if the difference of the total semileptonic decay rates between Dg 
and mesons could be measured experimentally. 

As will be discussed in Sec. 2 this contribution is described, as usual, by dimension 6 
four-quark operators. The matrix element of the relevant combination of the operators over 
the B mesons vanishes in the limit of factorization. However, also as usual, the considered 
term is greatly enhanced numerically, so that a violation of the factorization relation by 
10% makes this contribution about twice larger than the estimate of the corresponding 
uncertainty in Ref.0. Furthermore, the effect explicitly depends on the light quark flavor 
and thus is generally different in the decays of B^ and B^ mesons. In Sec. 3 it is argued 
that a deviation from factorization of such magnitude would not be unusual, and in fact is 
possibly indicated by the observed difference of the lifetimes of the Dg and mesons, if 
that difference is analysed within the same OPE framework. Moreover, the discussed effect 
in the constrained decay rate oi B ^ Xud-v is very closely related to a similar effect in 
semileptonic charm decays, in particular to the yet unmeasured difference of total inclusive 
semileptonic decay rates of Dg and D^. Thus a quantitative understanding of the discussed 
contribution would become possible, once this difference of the semileptonic decay rates is 
measured experimentally. The concluding Section 4 contains discussion and summary. 
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2 The third term in OPE for constrained inclusive de- 
cay rate. 

The optical theorem of the scattering theory relates the total decay rate Th of a hadron 
Hq containing a heavy quark Q to the imaginary part of the 'forward scattering amplitude'. 
For the case of inclusive decays generated by a particular term in the weak interaction 
Lagrangian the latter amplitude is described by the following effective operator 



Leff = 2Im 



(1) 



d''xe'P''T{Lly{x),Lw{0)} 
in terms of which operator (at = mg) the total decay rate is given by[| 

= {Hq\ Leff \Hq) . (2) 

Using in eq.(P the term 

Luh = —^yf^ 7m (1 - 75) b) (3) 

with = {i'jfj, (1 — 75) z/) in place of L\y, one would find the total inclusive decay rate 
of -B — s> Xuity. The effective operator (|l]) is evaluated using short- distance OPE. The 
leading term in the expansion describes the perturbative decay rate, while subsequent terms 
containing operators of higher dimension describe the nonperturbative contributions. In 
particular, the leading term with no QCD radiative corrections applied, gives for the total 
inclusive rate of the decays B i v the familiar expression 

^° - 192 773 ■ 

For a constrained inclusive rate in the case of semileptonic decays and with selection in 
the invariant mass of the lepton pair, one should rather consider, instead of the full correlator 
(2), a correlator of the quark currents, and then perform integration over the lepton pair 
phase space with the weight function corresponding to the considered selection of events [0. 
However for the third term in the OPE, suppressed by and expressed through four- 
quark operators of dimension 6, this procedure is not really necessary, since formally this 
contribution corresponds to an infinitesimally small momentum carried away by the final 
light quark, and all the momentum of the initial heavy quark flowing through the lepton 
pair. Thus this contribution formally comes only from = mg, and does not depend on the 
^ The non-relativistic normalization for the heavy quark states is used here: {Q\Q^ Q\Q) = 1- 
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lower cutoff q^. (This behavior of the discussed nonperturbative effect, first pointed ou in 
Ref. 1^, was also used in the estimate of Ref.p].) Therefore the effect of this term constitutes 
a fixed fraction of the total rate of the decays B — > Xu i v- 

It is quite straightforward to find the expression for the third term in the expansion of the 
correlator (|Tp appropriate for the decays B X„ iv |]^, in terms of four-quark 

operators normalized at = m;,: 

(3) _ 2G\\Vuh\^ml / X 



where the following notation is used for the relevant four-quark operators: 



0^_^ = {hLl^,qL){qLl^,hL) , 0^_p = {hRqL){qLhR) , 

n-A = iht''l,qL)iqLt''l,bL) , n_p = ibnt''qL)iqLt^bR) . (6) 

(The operators T, containing the color generators t", will appear in further discussion.) 

The matrix elements of the four-quark operators over the B mesons are parameterized 
g as 

{B\0-_^\B) = ^^B,, {B\0-s-p\B) = ^-^B,, (7) 

where is the meson annihilation constant, and Bi and B2 are generally phenomenological 
parameters, "bag constants". In the limit of naive factorization, i.e. where the product of 
bilinear operators is saturated by only the vacuum insertion, these parameters are Bi = 
B2 = 1 for the B meson containing the same light quark as the operator {B^ mesons for the 
case of the operators O"), and are zero for the other two mesons {B^ and Bg in this case). 
Clearly, in this limit the matrix element of the four-quark operator in eq. (|^) over either of the 
B mesons is vanishing^. However, as useful as the factorization relation might be for other 
estimates, one does not expect this relation to be valid with better than about 10% accuracy, 
and it is conventionally implied that deviations of such order of magnitude are present. (For 
numerical estimates of non-factorizable terms see e.g. |TB|, IT^, 0.) Moreover, it should 



be emphasized that these deviations refer not only to the relation between the constants: 
Bi = B2 for the meson with the same flavor as that of the quark q in the operator O'^, but 
that these constants are generally non zero when the light quark flavors do not match [^. 
In other words, there are separate sets of the constants B describing the matrix elements 



^The naive factorization for the discussed effect actually correspond to the annihilation decay iv, 
which obviously is forbidden by chirality in the limit where the lepton mass is neglected. 
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Bq/q = {Bqi\0''\Bqi) for q' = q and for q' ^ q. Taking into account also the violation of the 
flavor SU(3) would further proliferate the constants B. 

In any case, a violation of the relation Bi = B2 by about 10% would not be surprising at 
all. On the other hand, a deviation from the exact factorization of such moderate magnitude 
would be quite significant for the discussed effect in the B — ^ i v decays. Indeed, in terms 
of the "bag constants" the contribution bV to the decay rate of the term @ is expressed 
according to equations (H) and (0) as 

n2 IT/ 12 f2 ™3 

bV ^ ^^'^""l ^^""^ (52 - . (8) 
12 vr 

Comparing this with the dominant semileptonic decay rate of -B — Xc £ z/, and taking into 
account the known semileptonic branching ratio Bsi{B), one estimates that the discussed 
nonperturbative contribution corresponds in terms of the branching ratio to the correction 

SBiB - X„f « 3^9 [^)' (^) IK.P . (9) 

which easily can exceed by a substantial amount the previous estimate of this correction 
as 2.0 \ Vub\'^ if the deviation from factorization is, as expected, at a 10% level. 



3 Non-factorization effects in D mesons 

The deviation from naive factorization in mesons can in fact be probed more quantitatively 
from the data on charmed D mesons. Indeed, in the limit of factorization, the OPE for the 
dominant Cabibbo unsuppressed nonleptonic decays due to c ^ sud predicts equal rates 
of decay for Dg and mesons [^, 10 1. Experimentally it is known by now |l]] that there 



is a noticeable difference in lifetimes of these mesons: t{Ds)/t{D^) = 1.20 ± 0.025, which 
cannot be described by spectator dependent effects in Cabibbo suppressed decay channels, 
or by the flavor SU(3) breaking Although this discrepancy can be merely attributed 
to the overall inaccuracy of the OPE in the inverse of the charmed quark mas^, a more 
constructive approach would be to attempt describing this difference in lifetimes as due to 
deviations from factorization (see also in f^, |^). As will be discussed further in this section 
such approach also can be tested by measuring the difference of inclusive semileptonic decay 



•^In this respect the situation is no better for the expansion of a constrained inclusive rate of the decays 
B — > Xu £ where the expansion is governed by jic < nic M. 
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rates of the Dg and D mesons, which difference is somewhat more directly related to the 
discussed matrix elements over the mesons in equations (^) and (^. 

In the limit of flavor SU(3) symmetry the difference of the dominant inclusive nonleptonic 
decay rates of and Dg mesons can be written in terms of matrix elements of four-quark 
operators (normalized at /i = m^) as 

r(D°) - r(D,) = ^^^'^^crnlfDrriD _ ^ns _ 3 ^ 3 ^ ^^^^ 

where 6c is the Cabibbo angle, C+ and C_ are the well known short- distance QCD renormal- 
ization coefficients for nonleptonic weak interaction: C_ = C^'^ = {as{mc) / as{mw)y^^'^^ , 
and the flavor non-singlet coefficients B and e parameterize the following differences of the 
matrix elements: 

- {Oy-A " Ov-a)ds-do 

- {Og-p - Os_p)ds-D0 

/n~'S rpu \ 

- \J-V-A ~ J-V-AlDs-D^ 

where the operators O and T are the same as in eq.(^ with the h quark being replaced by 
c and the notation {X)a~b = (^1-^1^) ~ {B\X\B) is used. (The parameters ei and 62 both 
vanish in the limit of factorization.) It should be also mentioned that no attempt is being 
made here to allow for the breaking of the flavor SU(3) symmetry, thus no distinction is 
made between the annihilation constants or masses of the Dg and Dq mesons. 

The expression (|TU]) for the difference of the total decay rates corresponds numerically to 

r(i5°) - TiDg) ^ 3.3 J^^^ ) (i?r - Br -l^r + l ^r) ps-' ■ (12) 

Comparing this estimate with the experimental value for the difference of the total decay 
rates: 0.41 ± 0.05 one arrives at an estimate of corresponding combination of the non- 
singlet factorization parameters: 

- -^e'l' +^er ^0.12 , (13) 

which perfectly complies with the understanding that non-factorizable contributions are at 
a level of about 10%. 
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The estimate (|T3D of the non-factorizable terms however can serve only as a semi- 
quantitative indicator of the magnitude of the spectator effects in the inclusive rate of the 
processes B X^fi/ described by a different combination of the factorization parameters 
in eq.(^) than in eq. ([TB|) . A somewhat more direct test of the relevant combination of the 
parameters would be possible from the difference of the total semileptonic decay rates of Dg 
and mesons. Indeed, in the limit of flavor SU(3) symmetry this difference arises only in 
the decays due to c ^ s iu and is given in terms of the operators normalized at /x = rric as 

Gl cos2 9c ml fl niD _ 

(14) 

Given that the total semileptonic decay rate of the meson is approximately 0.16ps~^ [0, 
the discussed difference can easily amount to a quite sizeable fraction of the semileptonic 
rate, provided that 5"'' — 32^^ ~ 0.1. A more precise calculation of the relative effect of 
this difference in terms of the parameters B'^^ and Bl^" is somewhat difficult to achieve 
at present. The problem here is that the observed semileptonic decay rate of already 
includes, in comparison with the 'bare' parton rate, quite substantial negative QCD radiative 
corrections as well as the 0(m~^) corrections |T^, while neither of these is included in eq . ([T^) . 
If compared with the 'bare' parton semileptonic decay rate of c — s £ z/, the discussed effect 
is of a somewhat more moderate relative magnitude: 

Vsl{D^) -Tsl{.Ds) _ ^ , ( fp Y tuns j^ns^ p-N 

In either case however the difference should be quite conspicuous and can amount to several 
tens percent, provided that — i^g''! ~ 0.1. 

A measurement of the difference of the inclusive semileptonic decay rates of the and 
Dg mesons would make it possible to more reliably predict the difference of the corresponding 
decay rates between B^ and mesons: T{B^ — > Xu i v) —V{B^ — > i v\ which, according 
to the previous discussion, is dominantly concentrated in the upper part of the spectrum 
of the invariant mass of the lepton pair 0, §• At the level of accuracy of the present 
discussion the only difference between the theoretical expressions for B and for D mesons 
arises through a different normalization point of the four-quark operators in the equations 
(H) and (0). Taking into account the 'hybrid' evolution of the operators containing h quark 
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down to /i = rric gives the relation between the non-singlet factorization constants: 

Brim,) - B^im,) = ^^^-^ [Brim,) - BTim,] - ^ ^"'^ ~ [eTim,) ' eTim,] , 

(16) 

where k = (as('^c)/«s('^fe))- However, modulo the unlikely case that the difference of the 
constants e in this relation is much bigger than the difference between the constants B, the 
renormalization effect is quite small, and most likely is at the level of other uncertainties in 
the considered approach (such as the accuracy of the flavor SU(3) symmetry, higher QCD 
corrections, contribution of higher terms in m~^, etc.). Thus with certain reservations, one 
can use the approximate relation Br{mb) —B'r{iTib) ^ Br{mc) —BriiTic) to relate directly 
the differences in the inclusive semileptonic decay rates: 

r(B°^x„£z/)-r(i?±^x.£z/)^|K.lVlKsP^4 (r.K^°) - r.z(i^s)) • (ir) 



4 Discussion and summary 



Naturally, the largest uncertainty of the approach pursued in the present paper is perceived 
to be coming from applying the operator product expansion for the amplitude in eq.([l|) 
in Minkowski space at energy equal to the mass of the charmed quark. This raises both 
the issue of the contribution of higher terms as well as that of the applicability of quark- 
hadron duality at such energy. The situation is further complicated by the poor knowledge 
of the hadronic matrix elements of four-quark operators, including the deviation from the 
naive factorization in mesons. Same uncertainties are pertinent in the same, or greater, 
measure to a calculation of the constrained inclusive rate of the decays B X^iv at 
> [Mb — MoY- At present these issues are unlikely to be resolved purely theoretically, 
and an experimental input is acutely needed. The accuracy of the approach can be probed 
by testing at a quantitative level its theoretical predictions, most notably that of a large 
enhancement |TT| of the semileptonic decay rates of the strange charmed hyperons S^, S° as 



compared to that of Ac. It was argued previously [H, ^ that the non-factorizable terms give a 
sizeable contribution to the overall inclusive semileptonic decay rate of D mesons. However, 
given the uncertainty in other parameters (the charmed quark mass, higher QCD corrections, 
etc.) it could be somewhat ambiguous to assess the light-flavor singlet effect of these terms in 
the decay rates. As discussed in this paper, an experimental measurement of the difference 
of inclusive semileptonic decay rates between and D^, determined by the non-singlet part 
of the non-factorizable terms, is more likely to shed light on the problem of the deviations 
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from factorization in heavy mesons. Unless a better understanding of the non-factorizable 
terms is achieved, they will stand in the way of determining the mixing parameter 
from inclusive semileptonic decay spectra with accuracy better than 0(15%), including the 
method discussed in Refs.[^, 0. 

In summary. It is pointed out that non-factorizable terms can provide a quite noticeable 
contribution to inclusive rates of the decays B Xud-V-, which depends on the flavor of 
the light spectator quark, and is thus different for and B^ mesons. This yet unknown 
contribution can limit the theoretical accuracy of determining \ Vub\^ from the lepton spectra 
beyond the kinematical limit for decays with charm in the final state. The effect of such 
contribution can be substantially larger than the previously estimated uncertainty, if 
the non-factorizable terms are evaluated from the experimentally observed difference of the 
lifetimes of Ds and Dq mesons. A more direct conclusion about the non-factorization effects 
in charmless semileptonic decays of the B mesons can be drawn once the difference of the 
total semileptonic decay rates of D*^ and Dg is measured experimentally. If the current 
understanding is correct and the non-factorizable terms amount to about 10%, the latter 
difference can amount to a significant fraction of the total semileptonic decay rate of charmed 
mesons. 
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